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Abstract. Initial results from simulations of 4-fluid MHD turbulence in molecular 
clouds are presented. The species included in the simulations are ions, electrons, nega- 
tively charged dust grains and neutrals. The results indicate that, on length scales of a 
few tenths of a parsec, multifluid effects have a significant impact on the properties of 
the turbulence. In particular, the power spectra of the velocity and magnetic fields are 
significantly softened, while the PDF of the densities of the charged and neutral fluids 
are appreciably different. Indeed, the magnetic field strength displays much less spatial 
structure on all lengthscales up to 1 pc than in the ideal MHD case. The assumptions 
of ideal MHD therefore appear to be inadequate for simulating turbulence in molecular 
clouds at these length scales. 



1. Introduction 



It is g enerall y believed that mole cular clouds are turbulent (see the revie ws of M ac Low & Klessen 



2004; Elmeg reen & Scald l2004h . Observationally (e.g. lLarsonlll98lh this turbulence 
appears to be supersonic, with RMS Mach numbers of anything up to 10, and either 
super-Alfvenic or trans- Alfvenic. The amplitude of this turbulence makes it likely to 
be an important ingredient in our understanding both of th e dynamics of molecular 
clouds and of the process of star formation dElmegreenll 1993b and hence developing an 
understanding of this phenomenon is of considerable interest. 

Many authors have addressed the issue of MHD turbule nce in the context of 
molecular clouds using both the ideal MHD approximation (e .g. Mac Low et al. 19981: 
Mac Lowll999l : IOstriker et aLll200ll : lGiover & Mac Lowl2007l:lLemaster & StonefcOO 
and, more recently, various flavours of non-ideal MHD dOishi & Mac Lowl2006l : iLi et al] 
20081 : iKudoh & Basull2008l : iDownes & Q'Sullivanll2009ll201lh . Ambipolar diffusion 
has been found to cause greater temporal variability in the turbulence statistics. Clearly 
of lesser significance, the Hall effect is in principle capable of inducing topological 
changes in the magnetic field which are quite distinct to any influence caused by am- 
bipolar diffusion. 

The non-ideal MHD effects caused by multifluid physics occur on relatively small 
length scales (fractions of a parsec). However, since three dimensional MHD turbu- 
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lence involves the transfer of energy from the energy injection scale to ever smaller 
scales until the dissipation length scale of the system is reached, it is virtually inevitable 
that multifluid effects will have an impact on some part of the energy cascade. If the en- 
ergy transport is non-local in fc-space then these effects may influence the entire energy 

cascade. 

Downes & O'Sullivanl ( 200911201 lh performed simulations of decaying, non-ideal 



MHD and multifluid MHD, molecular cloud turbulence incorporating parallel resistiv- 
ity, the Hall effect and ambipolar diffusion. They found that the Hall effect has sur- 
prisingly little impact on the behaviour of the turbulence: it does not affect the energy 
decay at all and has very limited impact on the power spectra of any of the dynamical 
variables, with the exception of the magnetic field at very short lengthscales. 

In this paper we use the HYDRA code (lO ' Sullivan & DowneslboOd 120071) to in- 
vestigate driven, isothermal, multifluid MHD turbulence in a system consisting of 1 
neutral and 3 charged fluids. The aim of this work is to determine the influence of 
multifluid effects on the behaviour of driven turbulence in molecular clouds. An added 
benefit of our properly multifluid approach is that we can self-consistently deduce the 
behaviour of the charged spec ies and, in principle, make links with observations such 
as those of Li & Houdd (I2008h . 



2. The model 

Molecular clouds are, to a great extent, weakly ionised systems. This allows us to make 
several simplifying assumptions, as follows: the bulk flow velocity is the neutral veloc- 
ity; the majority of collisions experienced by each charged species occur with neutrals; 
the charged species' inertia is unimportant and, finally, the charged species' pressure 
gradient is unimportant. Ma king these assumptions allows us to derive a relationship a 
generalised Ohm's law (e. g. IFalldl2003l : Idolek & Robergdl2002l) . and this allows us to 
avoid having to calculate the electric field explicitly from the charge distribution. 

The resu lting weakly ionised, multifluid MHD equations are a nonlinear system 
of PDEs (see lCiolek & Robergdl200l : lFaUdl200l : lQ' Sullivan & Downesll2007h and, in 
general, must be tackled using simulations. The method used by HYDRA to integrate 
these equations has been elucidated previously in the literature (e.g. lO' Sullivan & Downes! 
I2006l . l2007h and we refer the reader to these papers for a more in-depth discussion. The 
method is explicit and scales extremely well on massively parallel systems, displaying 
strong scaling from 8192 up to 294,912 cores with 73% efficiency on the JUGENE 
BlueGene/P system. 



3. Numerical set-up 

We run our turbulence simulations in a cube of side 1 pc with periodic boundary con- 
ditions on all faces. The density and magnetic field are initially uniform with val- 
ues of 10 6 cm -3 and 1 mG respectively. The (isothermal) sound-speed is chosen to be 
3.3 x 10 4 cm s -1 or, equivalently, the temperature is 30 K. The four fluids in the system 
are initially at rest and the properties of the charged species are chosen to approximate 
those of dust grains, ions and electrons with the dust grains making up 1% of the mass 
density of the neutrals. We require charg e neutrality and that the resistivities in the 
generalised Ohm's Law be as described in IWardle & N d (Il999h . The grid used is uni- 
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Figure 1 . Evolution of the RMS Mach number with time (left panel) and the time- 
averaged compensated neutral velocity power spectra (right panel) for each simula- 
tion in the resolution study. 



form and made up of iV| grid points with N g points in each of the x, y and z directions. 
We emphasise here that this work is directed at understanding the nature of multifluid 
MHD turbulence. While a physical system with the properties just described may well 
be subject to dynamically important self-gravity, it seems reasonable that the optimal 
way of approaching understanding the behaviour of molecular clouds is to attempt to 
understand the turbulence itself before including the effect of self-gravity which will, 
undoubtedly, significantly impact the behaviour of the system. 

To drive turbulence we add velocity increments to the neutral veloc ity at each time- 
step. T he incremental velocity field, Su, is defined in a manner similar to lLemaster & Stone 
(120091) . Each component of 6u is generated from a set of waves with wave numbers, 
3 < k < 4 where k = |k|. The amplitudes of these waves are drawn from a Gaussian 
random distribution with mean 1.0 and deviation 0.33 while the phases of the waves are 
drawn from a uniform distribution between and 2n. A new realisation of the incre- 
mental velocity field is generated after 10~ 3 t c has elapsed since the previous realisation 
was generated. 

The rate of injection of energy is the same for all simulations in this work and is 
given by — jj-j = 300 where po is the initial (uniform) density, L is the length of side 

of the computational domain and c s is the sound speed. In our simulations this energy 
injection rate yields steady state turbulence with an RMS Mach number of around 5 and 
an Alfvenic Mach number (with respect to the unperturbed magnetic field) of around 
1.2. 

The nomenclature used here to refer to the simulations is defined as follows: each 
simulation is referred to as either mf-xxx or mhd-xxx with "mf" standing for multifluid 
and "mhd" standing for ideal MHD. The "xxx" gives the value of N g . 



4. Resolution study 

Figured] contains the Mach number as a function of time and the compensated power 
spectra of the neutral velocity for three simulations (mf-128, mf-256 and mf-512) in our 
resolution study. While there are differences between mf-128 and mf-256 it is clear that 
the differences between mf-256 and mf-512 are minor, indicating that both simulations 
resolve the multifluid physics occurring in this system reasonably well. 
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Figure 2. Plot of the RMS Mach number (left panel) and the perturbed magnetic 
energy (right panel) as a function of time for mf-512 and mhd-512. 



5. Results 



In this section we analyse some aspects of the turbulence generated in our simulated 
systems. We compare the properties of the mf-512 simulation with that of mhd-512 
in order to gain a sense of the differences introduced by multifluid effects at the 1 pc 
level. The left panel of Figure Q] indicates that for t > 0.12 t c the turbulence has reached 
a quasi-steady steady state. Any time-averaged quantities discussed in this section are 
averaged over times after t = 0.12 t c . 



5.1. Energy dissipation 



Figure El contains plots of the behaviour of the mean, mass-we ighted RMS Mach num- 
ber (see Lem aster & Stonell2009l : lDownes & Q'Sullivanll201ll) and the perturbed mag 



netic energy as functions of time. It is clear that the Mach number reached in each case 
is very similar and in fact the difference in time-averaged RMS Mach number after 
t = 0.12 is less than 1%. Therefore, at these length scales and with this numerical res- 
olution, we do not detect any difference in the energy dissipation rate due to multifluid 
effects. It is worth noting, however, that the perturbed magnetic energy in mf-512 is 
lower than that in mhd-512 (Figure [2l right panel) which indicates that the energy in 
the turbulence is less weighted toward the magnetic field in a multifluid MHD system 

than in an ideal MHD one. 

The work of iDownes & O' Sullivan! ( 20091 . l201lh implies that the energy dissipa- 



tion rate is higher in multifluid MHD systems. Their simulations were carried out in a 
cube of length 0.2 pc and the decay rate of the turbulence is undoubtedly higher in their 
multifluid simulation than in the ideal case. While this point is worthy of further study, 
it may be the case that the difference in the length scales, which determines the physical 
resolution of the simulations, may be partly responsible for the apparent di screpancy in 
the results presented here and those of lDownes & Q' Sulliv an ( 20091. l201ll) . 



5.2. Power spectra 

We now turn to the power spectra of the various quantities of interest in these sim- 
ulations. The nature of these spectra will give us some idea of the influence of the 
multifluid effects on the turbulence. We time-average all the power spectra discussed 
here over the time-interval [0.12 t c : 0.22 f c ]. 
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Figure 3. Time-averaged, compensated power spectra of the velocity field (left 
panel) and the magnetic field strength (right panel) for mhd-512 and mf-512. The 
velocity spectra are compensated with k l 5 while the magnetic field spectra are com- 
pensated with k 2 . 

Figure [3] contains plots of the velocity and magnetic field strength power spectra 
for mhd-512 and mf-512. It is clear that the velocity field in mhd-512 has less structure 
at length scales below about 0.05 pc than that of mf-512. The difference between the 
two systems is more dramatic in the power spectra of the magnetic field strength where 
simulation mf-512 displays considerably less structure at virtually all length scales be- 
low the driving scale. This is, perhaps, what we would expect given the fact that the 
multifluid effects manifest themselves most directly in the induction equation. 

The left panel of figure |4] contains plots of the power spectra of the densities for all 
of the fluids in mf-512. At these length scales the electrons and ions have almost identi- 
cal power spectra while the neutral fluid has considerably more power at all lengthscales 
below ~ 0. 1 pc. These differences are due to the fact that the electrons and ions are more 
closely tied to the dynamics of the magnetic field than are the neutrals. The magnetic 
field, as a result of ambipolar diffusion, has less power at shorter length scales in par- 
ticular (see Figure [3]), and hence the species closely tied to it also have less fine-scale 
structure. Interestingly the density PDFs of the fluids in mf-512 (right panel of figure 01) 
indicate that, though all fluids are distributed in an approximately lognormal way, the 
distributions of the charged species are considerably more peaked and narrower. This 
is another manifestation of the effect just noted: since the magnetic field is somewhat 
decoupled from the bulk (neutral) flow it is, in particular, not compressed or rarefied by 
the turbulent flow as much as the neutral flow. This means that the species tied to the 
magnetic field will not be compressed or rarefied as much as the neutrals and hence the 
PDF of their densities will have a lower standard deviation than that of the neutrals. 



6. Conclusions 

In this work we present 4-fluid MHD simulations of turbulence in molecular clouds, 
neglecting self-gravity. The aim is to investigate the impact of multifluid MHD effects 
on the turbulence itself while the inclusion of self-gravity will be the subject of future 
work. The results indicate that, on the lengthscales at which the simulations were run 
(~ 1 pc) multifluid effects have an appreciable effect on both the velocity and density 
power spectra, and there is a significant difference between the PDFs of the neutral and 
charged species. 
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Figure 4. Time-averaged, compensated power spectra of the densities of the var- 
ious fluids in mf-512 (left panel) and the probability distribution functions of the 
logarithm of their densities (right panel). 



This all implies that, even on fairly large length scales, multifluid effects are dy- 
namically important in molecular cloud turbulence. 
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